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Although the notion of entropy lies at the core of statistical mechanics, 
it is not often used in statistical mechanical models to characterize phase 
transitions, a role more usually played by quantities such as various order 
parameters, specific heats or susceptibilities. The relative entropy induces 
a metric, the so-called information or Fisher-Rao metric, on the space of 
parameters and the geometrical invariants of this metric carry information 
about the phase structure of the model. 

In various models the scalar curvature, TZ, of the information metric 
has been found to diverge at the phase transition point and a plausible 
scaling relation postulated. For spin models the necessity of calculating in 
non-zero field has limited analytic consideration to one-dimensional, mean- 
field and Bethe lattice Ising models. We report on previous papers in which 
we extended the list somewhat in the current note by considering the one- 
dimensional Potts model, the two-dimensional Ising model coupled to two- 
dimensional quantum gravity and the i/iree-dimensional spherical model. 
We note that similar ideas have been applied to elucidate possible critical 
behaviour in families of black hole solutions in four space-time dimensions. 

PACS numbers: 02.50.-r, 05.70.Fh 
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0. Introduction 

The notion of entropy, 

S = -^Paln{pa) , (1) 

a 

lies at the heart of statistical mechanics where pa is the probability of being 
in state a. It gives a measure of the number of microstates accessible to a 
system = InJl. It is also possible to define a relative entropy with respect 
to some reference set of configurations of probability (such as those at 
the critical point or zero temperature), 

G{p\r) = -J2p^ln(P^). (2) 

If 9 represents parameters characterizing the class of models under consider- 
ation, the relative entropy, G, induces a metric for two "close" configurations 
e ande + 59: 

= G {p{9)\p{9 + 59)) . (3) 
Since G{p\p) = and dG/d9 = we find 

d^G 

% 

If we assume a Gibbs-type distribution for the p, 



P,. = |exp(-9,$i) ^S= +lnZ, (5) 

then the metric may be written as 

This metric is well known in statistics as the Fisher-Rao metric, and charac- 
terizes the "closeness" of probability distributions For a one-parameter 
distribution the expectation value is the Fisher Information, I [9), given by 



and is inversely related to the variance I(^)Var(^) = 1. In the multiparam- 
eter case this generalizes in the obvious manner to the covariance matrix. 
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Various authors |21IS1I3 have observed that the geometric invariants of 
the metric, in particular the curvature, TZ, might be used to characterize 
the phase structure of statistical mechanical models. The first observation 
to note is that a non-interacting model, the ideal gas, has a flat geometry 
since T^ideai = 0. The suggestion was then that the curvature was a measure 
of the strength of interaction, an observation confirmed by a calculation in 
the case of a van der Waals gas (which is interacting): 



4 a/3 F(a,/3) 

Q~r n^., m2' 



RvdW — g~7v ^ 

3 V D[a,p) 



where /3 is the inverse temperature, a the pressure, and v{a, (3) and -F(a, (3) 
are expressions appearing in the equation of state. This shows that the 
curvature diverges along the spinodal line where D{a,(3) = 0. 

These examples suggest that phase transitions are manifested as diver- 
gences in the scalar curvature associated with the Fisher-Rao metric, at 
least in two-parameter models. A natural question to ask is whether any 
of the standard scaling exponents associated with a transition can then be 
extracted from the behaviour of the curvature. This is best pursued by 
consideration of some specific examples. If we consider a spin model in 
field, the manifold of parameters, is two-dimensional and parametrised 
by {6^ , 9'^) = (/?, h). In this case, the components of the Fisher-Rao metric 
take the particularly simple form 

Gi, = d^djf, (9) 

where / is the reduced free energy per site and di = d/d9^. With the metric 
given in equ. ©, TZ may be calculated succinctly as 



7^ 



2G2 



d}f dpdnf dlf 

djf djd,f dpdlf 

dldnf dpdlf dlf 



(10) 



where G = det(Gjj). Since the only scale present near criticality for a 
model displaying a higher-order transition is the correlation length, ^, it 
has been hypothesised on dimensional grounds that TZ ~ where d is the 
dimensionality of the system (2* J, 4 . If we assume that hyperscaling holds, 
ud = 2 — a, this leads to 

|^|(2-a)/^_ (11) 

To test the behaviour of 7^, one requires models which are solvable in 
field in order to obtain analytic expressions, and these are rather thin on 
the ground. Indeed, TZ has been calculated for the mean-field and Bethe- 
lattice Ising models [S] and the above scaling behaviour verified. It has also 
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been calculated for the one-dimensional Ising model where it takes the 
remarkably simple form, 

^ , cosh/i 

Rising = 1+ , . (12) 

V smh'' h + e-^^f^ 

In this case IZ is positive definite and diverges only at the zero-temperature, 
zero-field "critical point" of the model. The correlation length is given by 

= -ln(tanh(^)), (13) 

so that ^ ~ exp(2/3) near criticality, and (jllj) holds there with a = 1, = 1 
as expected^. 

Given this rather short list of examples, any further additions would 
be very worthwhile in order to see which features are generic and which 
are particular to the models concerned. We expand the list incremen- 
tally here, looking at the one-dimensional Potts model, the too-dimensional 
Ising model coupled to two-dimensional quantum gravity and the three- 
dimensional spherical model, which shares the critical exponents of the just 
mentioned coupled Ising model. We also discuss the application of similar 
ideas to critical behaviour in various families of black holes. The present 
paper is essentially an amalgam of results presented in pi [71 IH] , with some 
added discussion of related work in black hole physics. 

1. One 

The partition function for the ID (^-state Potts model is given by 



ZN{y,z) = ^exp 



N . 1\ ^ / 1 

13 J2 Si^j ,a,+,)--)+hJ2(6{a,,l)-- 



1^ .=1 



(14) 

where the spins, Uj G {1, . . . , g}, are defined on the sites, j G {1, . . . A^}, of 
the lattice and where we have defined y = exp(/3) and z = exp(/i) for later 
calculational convenience. The model may be solved by transfer matrix 
methods 0, just as the ID Ising model. For general q the full transfer 
matrix T{y, z) may be written as q — 2 diagonal elements, (y — l)(?/z)~^/'', 
and a 2 X 2 factor t{y, z): 



^ The Bethe lattice model also satisfies the postulated scaling, although there are some 
subtleties coming from the exponent a being zero 
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The partition function is Z]\i{y,z) = TT:T{y,z)^ and the eigenvalues of 
T{y, z) are Xo,Xi,. . . , Aq_i, where 



Ao 1_ 1 
Ai / "2 



y{l+z)+q-2± J{y{l -z)+q-2Y + {q- l)4z 



{yz) 1 
(16) 

and A2, • • • \-i = {y — l)(y-z) ^^'^ ■ The reduced free energy per site in the 
thermodynamic limit, — > 00, is thus given by / = — In Aq- 

It is straightforward to use this expression for the free energy in equ. ()1U() 
to obtain the curvature, TZ. In the current notation we re-derive the expres- 
sion for the Ising model ^ as 

Tz ■ 1 + y(i + -2) ^^-j 

i/y^ — 2y'^z + y'^z'^ + 4z 
The expression for general q is similar in form to this, and is 
t:, At \ , B{q,y,z) 

7^Potts = A{q, y,z) + — , (18) 

where the coefhcients A{q, y, z) and B{q, y, z) are smooth functions of y and 
z and do not diverge for finite (physical) temperature or field. Furthermore 

r/(Q, y, z) = [y(l -z)+q-2f + {q- l)4z. (19) 

The expressions for A{q,y,z) and B(q,y,z) are very lengthy for general q 
(although still easily obtained) so we do not reproduce them here. 

In zero-field (z = 1) the expression for TZ is much more compact and is 
written for general q as 

^ ^ jy + q- + (g - 2)y - (g - 2)(g - 1)) 

(g-l)(2y + g-2)2 " ^ ' 

We see that as y ranges from 1 to cxd, TZ ranges from (4 — q)/{q — 1) to 00. 
In particular, the sign of the y = 1 (/? = 0) limit of TZ changes at q = 4, 
although the general morphology of 7^ as a function of y and z remains the 
same for all g > 2 as we shall see below. 

The correlation length for the one-dimensional Potts model is defined in 
a similar manner to that of the Ising model, 

r' = -ln(^], (21) 



^ There is a factor of two difference in the definitions of /3, h between the Ising and 
Potts notations coming from the different spin definitions. 
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SO ^ ~ ?/ for z = l,y — > oo. We thus retrieve the scahng of TZ for the one- 
dimensional Potts model expected from equ. namely 7^ ~ y as y — > (X). 
The exponents, as for the Ising model, are a = 1,1^ = 1. 

The general features of TZ at non-zero temperature and field are perhaps 
easiest seen in a contour plot as a function of y and z. In Fig. 1 we show the 
Ising (q = 2) case which has certain non-generic features. The it/i symmetry 
of the Ising model appears as a z — > l/z symmetry in the plot of TZ. In 
addition, one can see that TZ is positive for all y and z. The maximum of 
TZ for a given y value lies along the zero field line at z = 1. In Fig. 2, 7?. is 
plotted for the 3-state Potts model. We see that there is no longer a 2; — > 1/z 
symmetry and that TZ is no longer positive definite. 

Some years ago Lee and Yang JJ], addressed the question of how the 
singularities associated with field-driven phase transitions in Ising-like spin 
models on lattices arose in the thermodynamic limit. This was later ex- 
tended by various authors to other models and to temperature-driven tran- 
sitions I12j . Lee and Yang observed that the zeroes of the partition 
function for a spin model in a complex external field on a finite lattice 
would give rise to singularities in the free energy. In the thermodynamic 
limit these complex zeroes move in to pinch the real axis, signalling the 
onset of a physical phase transition. Typically, the loci of zeroes are lines in 
the complex field or temperature plane and when the endpoints of such lines 
occur at non-physical (i.e. complex) external field values they can be con- 
sidered as ordinary critical points with an associated edge critical exponent, 
usually dubbed the Lee- Yang edge exponent [TT] . 

The Lee- Yang zeroes for the one-dimensional Potts model on a periodic 
chain with N sites are given by the solutions ^ °f 

= (Ai)^ + (Ao)^ = ^ Ai = exp(^^)Ao (22) 

where Ao,i are the eigenvalues given in equ. ()16() and n is odd. In the 
thermodynamic limit the locus of zeroes is determined by |Ao| = |Ai| or 

i]{q, y, z) = [y(l -z) + q-2f + {q- l)Az = 0, (23) 

which can be satisfied for complex (in the q = 2 Ising case, purely imaginary) 
values of h. 

From these considerations, it is clear that TZ will also diverge as the locus 
of zeroes is approached for both Ising and Potts models if we allow ourselves 
the liberty of an analytic continuation of the field to complex h values once 
TZ has been calculated, since TZ = A + B/ y/fj and ^, B are finite as — > 0. 
The presence of the square root means that the divergence is characterised 
by an exponent a = —1/2 which is the Lee- Yang edge exponent for the 
one-dimensional Potts (and Ising) model 
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The status of these observations is a httle unclear to us, since the cal- 
culation of TZ has assumed a real metric geometry throughout and such an 
arbitrary continuation in the final expression might be rather dangerous. It 
is nonetheless interesting that the Lee- Yang edge transition is still visible 
as a divergence in TZ. 



2. Two 

The solution of the Ising model on an ensemble of (4-regular) or 
$^ (3-regular) planar random graphs was first presented by Boulatov and 
Kazakov (14j . who were motivated by string-theoretic considerations, since 
the continuum limit of the theory represents matter coupled to 2D quantum 
gravity. They considered the partition function for the Ising model on a 
single n vertex planar graph with connectivity matrix G^j: 

ZsingUG^ /3, /i) = ^ exp I /3 ^ G^^a^aj + /i ^ , (24) 

then summed it over all n vertex graphs {G"} resulting in 

Zn=Yl Z^^^x,{G'\l3,h) , (25) 

before finally forming the grand-canonical sum over graphs with different 
numbers n of vertices 

^-|(^^)"-». (^«) 

where c = exp(— 2/3). This last expression could be calculated exactly as 
matrix integral over N x N Hermitian matrices, 



Z = - In J V(j)i P(/>2 exp (^-Tx 



-c</.i02-|(eVf + e-V!) 



(27) 



where the N ^ oo limit is to be taken to pick out the planar diagrams and 
the potential appropriate for <I>'^ (4-regular) random graphs has been shown. 

When the matrix integral is carried out the solution is given in para- 
metric form by 

1 . /z\ 1 dt 1 dt , 2 



Z = -\n[- -- g(t) + — -g{tY, (28) 

2 \gj g Jo t 2g'^ Jo t 
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where the function g{z) is 

" ' 0^ + .^] (29) 



giz) = \^z^ + 



(1-Z)2 (1-^2)2 



and B = 2[cosh(/i) - 1]. 

In the thermodynamic hmit the free energy per site is given by 

where zq = zo{P,h) is the appropriate low- or high-temperature solution of 
g'{z) = 0. When h = this may be solved in closed form, and although the 
solution is not available explicitly for non-zero h it can still be developed 
as a power series in h around the zero-field solutions in order to obtain 
expressions for quantities such as the energy, specific heat, magnetization 
and susceptibility. It was found that the critical exponents were given by 
a = — l,/?=l/2,7 = 2, so the transition was third order with, intriguingly, 
the same exponents as the 3D spherical model on a regular lattice |15| . 

If we carry out a perturbative expansion for the high-temperature solu- 
tion, which is symmetric in h and hence a series in even powers, we find 

1 (u- l)(2n2 -2u-M) 2 

^0 — 1 77] TTI n 

{u - l)(2n2 -2u + l){iu^ - + lOu^ - 5u^ + 5u + 1) ^4 
+ 24(2n - 1)9 ^ '^^^^ 

where the coefficients in the series are most naturally expressed in terms of 
u = exp(— /?) = ^/c, as above. 

The expected scaling form of the various components of TZ for a generic 
spin model in field is discussed at some length in ^16, , and we now recapitu- 
late these results briefly for comparison with the specific results for the Ising 
model on planar random graphs. The starting point is the scaling form of 
the free energy per spin near the critical point, 

fie,h) = X~'f{eX^%hX-'^) , (32) 

where e = (5c — 13 and a^^ah are the scaling dimensions for the energy and 
spin operators. For e > 0, i.e., in the unbroken high-temperature phase, we 
can use standard scaling assumptions to write this as 

f{e,h) = e^/'^^i:+{he~->^'''^), (33) 

where '0+ is a scaling function and we also define A = l/a^^ and C = —ah/a^ 
for later convenience. In terms of the standard critical exponents A = 2 — a 
and A + C = j3. 
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This generic scaling form can now be substituted into equ. to find 
the scahng behaviour of the scalar curvature H10|) near criticality (i.e. h = 0, 
e^O), 

^ _ {A + 2C)[iA + 2C)-iA-2)] _^ 
2A{A-1)M0) 
or, translating back to the standard critical exponents, 

^-2(2-a)(l-a)V'+(0)' " ^^^^ 

The discussion in ^Hl was intended to be as general as possible, one should 
note that for Ising-like models with a it/i symmetry all odd derivatives of 
the scaling function w.r.t. h will vanish so df^f = rather than e^~^^'~"ip_^_ (0). 
This does not affect the stated scaling relations. 

However, one feature of these scaling relations does have an impact 
on the observed scaling for the Ising model. Generically one expects that 
dpf = A{A — l)e^~'^ip^{0), which contributes to both the metric and the 
determinant involved in calculating TZ. If A > 2, i.e. a < 0, this naively 
suggests a vanishing d'^f at criticality, which will in general not be the case. 
There would instead be a contribution from a regular term, which would 
give a constant at the critical point. Having such a constant term, which 
we take to be A(A — l)(p(0) for notational convenience, modifies the scaling 
form of TZ in the case a < 0, ^ > 2 to 

(A + 2Cy , , , 

7^ = -A——^.^e-^ (36) 
2A{A - 1)0(0) ^ ^ 

or ^ 

^ " 2(2-a)(l-a),/.(0)' ' ' ^^'^^ 

so the critical exponent a no longer appears in the scaling exponent. 

By virtue of the Boulatov and Kazakov solution of the Ising model on 
planar random graphs [HI we can explicitly confirm these observations. 
Since a = —1,(3 = 1/2,7 = 2, we have A = 3, C = —5/2 and the modified 
discussion of scaling should apply. We can now take the series expansion 
for zq from equ. pij) . insert this into g{z) and use the resulting expression 
for / in equ. H30|) to calculate the various terms that appear in the scalar 
curvature TZ as power series in h"^. We find that the leading terms at /i = 0, 
with eu = u — Ucr = e/2 + . . . and Ucr = 1/2, and using /3, h as co-ordinates 
are 

'^- 704^- 

A glance back at equ. (|57j) shows that the modified scaling for A > 2 that 
these incorporate is, indeed, followed. 
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3. Three 

Berlin and Kac flT^ introduced the spherical model (and the Gaussian 
model) in an attempt to understand how generic some of the features of 
Onsager's solution iSj of the 2D Ising model are for ferromagnetic spin 
models, particularly for other dimensions. In the spherical model, the ±1 
condition on the value of the Ising spins is relaxed, whilst retaining a global 
constraint on the total spin magnitude. With Si denoting the value of a spin 
at a site z of a hypercubic lattice, the partition function is 

Z = j dsi...dsN exp l^YlsiSj + hY,s}j^(Y. " ^) ' (39) 

where is the total number of sites. This can be evaluated by exponen- 
tiating the constraint and using steepest descent, resulting in the following 
expression for the reduced free energy per site in the thermodynamic limit, 

N ^ oo: 

/ = ^ In + /3z - \g{z) + -J--r,. (40) 



4/3(z - d) 



where 



giz) = -j^;^ duJi... duJd In \^ - cos(wfc) j . (41) 

The saddle-point value of z which appears in the expression for the free 
energy in equ. (|40|) is determined from 

The solution reveals no transition for d = 1 and 2, and a transition with the 
exponents a = —1, /? = 1/2,7 = 2 for d = 3. 

It is useful to note that, with /i = 0, equ. gives 

dz 2 

and hence 

= (44) 

The critical point is given hy z = d = 3 and h = and the behaviour of 
g(z) in this region is determined by differentiating equ. (|41() twice and then 
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expanding for the small LVt values which give the dominant contribution. 
One finds 

A further differentiation gives 

g^'\z)^^{z-^)-'l^ (46) 
and an integration yields 

= -!-(.- 3)1/2 + (3), (47) 

where ^'(3) = (18 + 12^/2 - 10\/3 - 7\/6)[K(2^/3 + ^/Q - 2^2 - 3)]^ « 
0.505 462 019. . . is the massless 3D lattice propagator at the origin, which 
can be expressed in terms of one of the three classical Watson integrals 
jl9j and hence is given by the standard elliptic integral K{k'^). This latter 
expression can be combined with equ. ()42() with /i = to give 

(z-3) ~87r2(/3,-/?)2 (48) 

in which f3c = g'{2,)/2 k, 0.252 731009.... Equations ^ and ^ may 
then be substituted in equs. (|43|) and (|44jl to give the scaling of dz/dfi and 
(fz/dl3'^, 

dz 

lim — = lim{-4V27r(z - 3)^/2} = 0, 

which we shall employ below in the calculation of the scalar curvature. 

We now move on to examine the scaling of the various terms contributing 
to TZ in equ. (|10|) for the spherical model. As we have remarked, the h ^ —h 
symmetry in the free energy per site, /, of the spherical model means that 
any terms with an odd number of h derivatives will automatically be zero 
when h = 0, hence ffsh = fpph = fhhh = 0. This leaves the non-zero terms 
in the scaling region 
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c 



We see that the expected general scahng of each term (for a < 0) does 
indeed apply and that overall we have, as in equ. (|57|l. 

7e~e-2. (51) 

We thus see that calculating the scaling of TZ for the 3D spherical model for 
which a = —1 gives results in accordance with expectations from general 
scaling arguments which take into account the negative a, similarly to the 
Ising model on planar random graphs. 

4. Four 

The thermodynamics of black holes has been a subject of abiding in- 
terest since the pioneering work of Hawking 120] and similar ideas to those 
discussed in the previous sections for statistical mechanical models have also 
been applied to investigations of the critical behaviour of various families 
of black hole solutions in general relativity. Critical behaviour has arisen in 
several contexts in the study of black holes, ranging from the Hawking-Page 
j21j phase transition in hot Anti-de-Sitter space and the pioneering work by 
Davies on the thermodynamics of Kerr-Newman black holes, to the idea 
that the extremal limit of various black hole families might be regarded as 
a bona-fide critical point E31 [21 Hi iS] • 

It is the latter that is perhaps the closest to the work described here. In 
these studies the metrics used are the Ruppeiner metric 

gij = -didjS{E,N''), (52) 

and the Weinhold [27j metric 

gij = didjE{S,N''), (53) 

where 5 is the entropy and E the energy, with the A^" being other extensive 
variables such as scalar charges. Both metrics are related to the Fisher-Rao 
metric we have used by Legendre transforms of the appropriate variables, 
and have also been employed in a statistical mechanical context. 

For example , a Reissner- Nordstrom black hole has a Weinhold met- 
ric of the form 

d^w = ( ( 1 - ^ ) dS^ - SQdQdS + SSdqA , (54) 
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where Q is the charge, and we can see that one of the metric components 
vanishes at S = 2>Q^. This led Davies [22 to suggest that there was a 
phase transition at this point, but transforming to the Ruppeiner metric 
and choosing new co-ordinates gives a flat metric with i? = everywhere. 

On the other hand Kerr black holes which possess a spin J display a 
curved Ruppeiner geometry with a scalar curvature 



1 ^1-^-2 



that diverges at the extremal limit 



J 



M2 - ^''^ 

It would thus appear that the general framework of looking for the signal 
of phase transitions in the geometric invariants of the information metric 
(or related metrics) can also be employed in these circumstances. Indeed, 
the AdS / CFT correspondence identifies Hawking-Page type transitions with 
deconfinement in the dual gauge theories |24[ I28j so it might be profitable 
to explore the information geometry view of such transitions further. 
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Fig. 1. A plot of TZ for the one- dimensional Ising model for y = 1 ... 10, z = ... 5. 
The positivity of TZ and the expected z ^ 1/z symmetry are both apparent. 
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Fig. 2. A plot of TZ for the one- dimensional 3-state Potts model for j/ = 1 ... 10, 
z = . . .5. 7^ is no longer positive definite for physical values of y,z and there is 
no z ^ 1/z symmetry. In addition the maximum of TZ does not lie at 2; = 1. 



